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Abstract 

We present a simple recipe to construct exactly and quasi-exactly solvable Hamilto- 
nians in one-dimensional 'discrete' quantum mechanics, in which the Schrodinger equa- 
tion is a difference equation. It reproduces all the known ones whose eigenfunctions 
consist of the Askey scheme of hyper geometric orthogonal polynomials of a continuous 
or a discrete variable. The recipe also predicts several new ones. An essential role 
is played by the sinusoidal coordinate, which generates the closure relation and the 
Askey- Wilson algebra together with the Hamiltonian. The relationship between the 
closure relation and the Askey- Wilson algebra is clarified. 

1 Introduction 

For one dimensional quantum mechanical systems, two sufficient conditions for exact solv- 
ability are known. The first is the shape invariance p], which guarantees exact solvabil- 
ity in the Schodinger picture. The whole set of energy eigenvalues and the corresponding 
eigenfunctions can be obtained explicitly through shape invariance combined with Crum's 
theorem [2], or the factorisation method [3] or the supersymmetric quantum mechanics [I]. 



The second is the closure relation [5j. It allows to construct the exact Heisenberg operator 
solution of the sinusoidal coordinate r)(x), which generates the closure relation together with 
the Hamiltonian. The positive/negative energy parts of the Heisenberg operator solution 
give the annihilation/ creation operators, in terms of which every eigenstate can be built up 
algebraically starting from the groundstate. Thus exact solvability in the Heisenberg picture 
is realised. 

It is interesting to note that these two sufficient conditions apply equally well in the 
'discrete' quantum mechanics (QM) [6j [7J [8], [9] , which is a simple extension or deformation 
of QM. In discrete QM the dynamical variables are, as in the ordinary QM, the coordinate 
x and the conjugate momentum p, which is realised as p = —id x . The Hamiltonian contains 
the momentum operator in exponentiated forms e ±/3p , which acts on wavefunctions as finite 
shift operators, either in the pure imaginary directions or the real directions. Thus the 
Schrodinger equation in discrete QM is a difference equation instead of differential in ordinary 
QM. Various examples of exactly solvable discrete quantum mechanics are known for both of 
the two types of shifts [6j [7J [101 El [8j [9], and the eigenfunctions consist of the Askey-scheme 
of the hypergeometric orthogonal polynomials [HI [121 EE] of a continuous (pure imaginary 
shifts) and a discrete (real shifts) variable. 

It should be stressed, however, that these two sufficient conditions do not tell how to 
build exactly solvable models. In this paper we present a simple theory of constructing 
exactly solvable Hamiltonians in discrete QM. It covers all the known examples of exactly 
solvable discrete QM with both pure imaginary and real shifts [H [9] and it predicts several 
new ones to be explored in a subsequent publication |l4j . Moreover, the theory is general 
enough to generate quasi- exactly solvable Hamiltonians in the same manner. The quasi- 
exact solvability means, in contrast to the exact solvability, that only a finite number of 
energy eigenvalues and the corresponding eigenfunctions can be obtained exactly |15j . This 
unified theory also incorporates the known examples of quasi-exact ly solvable Hamiltonians 
[TBI H7J- A new type of quasi-exactly solvable Hamiltonians is constructed in this paper and 
its explicit examples will be surveyed in a subsequent publication [14] . One of the merits 
of the present approach is that it reveals the common structure underlying the exactly and 
quasi-exactly solvable theories. In ordinary QM, the corresponding theory was already given 
in the Appendix A of [5], although it does not cover the quasi-exact solvability. 

The present paper is organised as follows. In section two the general setting of the 
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discrete quantum mechanics is briefly reviewed and in §2.11 the Hamiltonians for the pure 
imaginary shifts and for the real shifts cases are given and the general strategy of working in 
the vector space of polynomials in the sinusoidal coordinate is explained. In §2.2[ based on 
a few postulates, various properties of the sinusoidal coordinate r)(x), which is the essential 
ingredient of the present theory, are presented in some detail. The main result of the paper, 
the unified form of the exactly and quasi-exactly solvable 'Hamiltonians', is given in §2.31 The 
action of the Hamiltonian on the polynomials of the sinusoidal coordinate is explained in 
§2.41 It simply maps a degree n polynomial into a degree n + L — 2 polynomial. Here L is the 
degree of a certain polynomial constituting the potential function in the Hamiltonian. The 
exactly solvable case (L = 2) is discussed in section three. In §3.11 the closure relation, which 
used to be verified for each given Hamiltonian, is shown to be satisfied once and for all by the 
proposed exactly solvable Hamiltonian. The nature of the dual closure relation, which plays 
an important role in the theory of discrete QM with real shifts and the corresponding theory 
of orthogonal polynomials of a discrete variable, is examined and compared with that of the 
closure relation in §3.21 The relationship between the closure plus dual closure relations and 
the Askey- Wilson algebra [HI [191 12Q1 [21] is elucidated in §3.31 In §3.41 shape invariance is 
explained and shown to be satisfied for the pure imaginary shifts case §3.4. II and for the real 
shifts case §3.4.21 The quasi-exactly solvable 'Hamiltonians' are discussed in section four. 
The QES case with L = 3 is achieved in §4.11 by adjusting the compensation term which is 
linear in r)(x). A new type of QES with L = 4 is introduced in §4.2[ which has quadratic in 
r)(x) compensation terms. It is shown that QES is not possible for L > 5 in §4.31 The issue 
of returning from the 'Hamiltonian' in the polynomial space to the original Hamiltonian Ti, 
is discussed in section five. This is related to the properties of the (pseudo-)groundstate 
<Pq. The final section is for a summary, containing the simple recipe to construct exactly 
and quasi-exactly solvable Hamiltonians. Appendix [A] provides the explicit forms of the 
sinusoidal coordinates with which the actual exactly and quasi-exactly solvable Hamiltonians 
are constructed. There are eight different tj(x) for the continuous variable x and five for the 
discrete x. Appendix [B] gives the proof of the hermiticity of the Hamiltonian, which is 
slightly more involved than in the ordinary QM. Appendix O recapitulates the elementary 
formulas for the eigenvalues and eigenvectors of an upper-triangular matrix, to which the 
exactly solvable (L = 2) 'Hamiltonian' in the polynomial space reduces. 
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2 'Discrete' Quantum Mechanics 



Throughout this paper we consider 'discrete' quantum mechanics of one degree of free- 
dom. Discrete quantum mechanics is a generalisation of quantum mechanics in which 
the Schodinger equation is a difference equation instead of differential in ordinary QM 
0, 0, HOI El E] . In other words, the Hamiltonian contains the momentum operator p = —id x 
in exponentiated forms e ±/3p which work as shift operators on the wavefunction 

e ±l3p i;{x) =i;{xTi(3). (2.1) 

According to the two choices of the parameter (3, either real or pure imaginary, we have 
two types of discrete QM; with (i) pure imaginary shifts, or (ii) real shifts, respectively. In 
the case of pure imaginary shifts, ij)(x =F ij), 7 £ K-^Oj we require the wavefunction to be 
an analytic function of x with its domain including the real axis or a part of it on which 
the dynamical variable x is defined. For the real shifts case, the difference equation gives 
constraints on wavefunctions only on equally spaced lattice points. Then we choose, after 
proper rescaling, the variable x to be an integer, with the total number either finite (N + 1) 
or infinite. 

To sum up, the dynamical variable x of the one dimensional discrete quantum mechanics 
takes continuous or discrete values: 

imaginary shifts : 168, x G (xi,x 2 ), (2.2) 
real shifts : i6Z, x G [0, N] or [0, oo). (2.3) 

Here X\, x 2 may be finite, — oo or +oo. Correspondingly, the inner product of the wavefunc- 
tions has the following form: 

imaginary shifts : (f,g)= / f*(x)g(x)dx, (2.4) 

J XI 

N oo 

real shifts: (f , g) = ^ f(x)* g(x) or ^ f(x)*g(x), (2.5) 

x=0 x=0 

and the norm of f(x) is ||/|| = a/ (/, /). In the case of imaginary shifts, other functions 
appearing in the Hamiltonian need to be analytic in x within the same domain. Let us 
introduce the *-operation on an analytic function, * : / i— >■ /*. If f(x) = a n x n , a n G C, 

n 

then f*(x) = 

a* n x n , in which a* is the complex conjugation of a n . Obviously f**(x) = f(x) 

n 
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and f(x)* = f*(x*). If / is an analytic function, so is g(x) = f(x—a), a G C. The ^-operation 
on this analytic function is g*(x) = (f(x* — a))* = f*(x — a*). If a function satisfies /* = /, 
then it takes real values on the real line. The 'absolute value' of an analytic function to 
be used in this paper is defined by = \/ f(x)f*(x), which is again analytic and real 

non-negative on the real axis. Note that the *-operation is used in the inner product for 
the pure imaginary shifts case (12.41) so that the entire integrand is an analytic function, too. 
This is essential for the proof of hermiticity to be presented in Appendix Q3] 

In quantum mechanics, the eigenvalue problem of a given Hamiltonian is the central 
issue. In this paper, we will consider the Hamiltonians having finite or semi-infinite discrete 
energy levels only: 

= E(0) < 8(1) < 8(2) < ■ ■ ■ . (2.6) 

Here we have chosen the additive constant of the Hamiltonian so that the groundstate 
energy vanishes. In other words, the Hamiltonian is positive semi-definite. It is a well known 
theorem in linear algebra that any positive semi-definite hermitian matrix can be factorised 
as a product of a certain matrix, say A, and its hermitian conjugate . As we will see 
shortly, the Hamiltonians of discrete quantum mechanics have the same property, both with 
the imaginary and real shifts. 

2.1 Hamiltonian and Strategy 

The Hamiltonian of one dimensional discrete quantum mechanics has a simple form 

Corresponding to the imaginary /real shifts cases, the parameter (3, the potential functions 
V±(x) and a sign factor e are 

imaginary shifts : (3 = 7, e — 1, V + (x) — V(x), V-(x) = V*(x), 

real shifts: (3 = i, e = -1, V+(x) = B(x), V-(x) = D(x), (2.8) 

with 7 G M^o- The potential function B(x) and D(x) are positive and vanish at boundaries: 

B(x) > 0, D(x) > 0, D(0) = ; B(N) = for the finite case. (2.9) 

As mentioned above, e ±/3p are shift operators e ± ' 3p /(x) = f(x =F ij3), and the Schrodinger 
equation 

HM^) =S(n)(f) n (x), n = 0,1,2,..., (2.10) 
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is a difference equation. The hermiticity of the Hamiltonian is manifest for the real shifts 
case because the Hamiltonian is a real symmetric matrix. For the imaginary shifts case, see 
Appendix [B] 

This positive semi-definite Hamiltonian (12. 7p can be factorized: 

H = A*A. (2.11) 

Corresponding to the imaginary /real shifts cases, A and A* are 

^ = z(e^/ 2 v ^(^)-e- 7p/2 v / V^)), A ] = -i( ^/V{x) e< p/2 - v^Re"^ 2 ) , (2.12) 
A= y f B{x) -e a ^/D(xj, A ] = ^B{x) - ^D{x) e- a . (2.13) 

The groundstate wavefunction (f>o(x) is determined as a zero mode of A, 

A(p ( x ) = 0- (2.14) 

The similarity transformed Hamiltonian Ti, in terms of the groundstate wavefunction </> has 
a much simpler form than the original Hamiltonian 7i: 

H = ^(x)- 1 oHo(j) (x) (2.15) 
= s(y + {x){e^ - 1) + V-{x){e-te - 1)). (2.16) 

In the second equation we have used (12.141) . 

In the following we will take 7i instead of 7i as the starting point. That is, we reverse the 
argument and construct directly the 'Hamiltonian' 7i (12.161) based on a certain function r](x) 
to be called the sinusoidal coordinate. The necessary properties of the sinusoidal coordinate 
will be introduced in the next subsection §2.21 The general strategy is to construct the 
'Hamiltonian' 7i in such a way that it maps a polynomial in rj(x) into another: 7iV n C 
Vn+L-2 C Voo- Here V n [n G Z> ) is defined by 

V n = Span[l, 77(x),..., 770r) n ], V M = lim V n . (2.17) 

n—>oo 

The goal is achieved by choosing very special forms of V±(x) as given in (I2.30I) - (I2.31I) . that is 
V(x) and V*(x) or B(x) and D(x) are polynomials of degree L in the sinusoidal coordinate 
rj(x) and its shifts r](x^i(3) divided by special quadratic polynomials in them. This provides 
a unified theory of exactly solvable and quasi-exactly solvable discrete QM. Exactly solvable 
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QM are realised by choosing Ji in such a way (L = 2) that 7iV n C V n is satisfied for all n. 
Then the existence of an eigenfunction, or to be more precise, a degree n eigenpolynomial, 
of Tt is guaranteed for each integer n. On the other hand, quasi-exact solvability is attained 
by adjusting the parameters of TC in such a way (L = 3,4) that H'Vm Q Vm is realised for 
an integer M. Here Ti 1 is a modification of 7i by the addition of the compensation terms. 
Then the 'Hamiltonian' Ji! has an M + 1 = dim(VA/)-dimensional invariant space, providing 
M + 1 eigenpolynomials of 7i'. After obtaining such (quasi-) exactly solvable 'Hamiltonian' 
H (H'), we have to find the (quasi-)groundstate wavefunction (j) in order to return to the 
true Hamiltonian H (H') by (I2.15p . It should be noted that the existence of such a (quasi- 
)groundstate wavefunction is not guaranteed a priori since we have started with 7i instead of 
7i. In the case of quasi-exactly solvable QM, the positive semi-definiteness of the Hamiltonian 
(12. 6p is in general lost due to the inclusion of the compensation terms to T~C'. 

2.2 sinusoidal coordinate 

Motivated by the study in [5], El E], let us define a sinusoidal coordinate r/(x) as a real (or 
'real' analytic r/*(x) = r/(x) in the case of pure imaginary shifts) function of x satisfying the 
following symmetric shift-addition property: 

7](x - if3) + r,(x + i/3) = (2 + r^ ] )ri{x) + r m x . (2.18) 

Here and r^\ are real parameters and we assume > —4. These two, and r^\, 
are fundamental parameters appearing in both exactly and quasi-exactly solvable dynamical 
systems. For the exactly solvable systems, these two parameters also manifest themselves 
(13. 4p in the closure relation, another characterisation of exact solvability, to be discussed in 
§3.11 Since a polynomial in r/(x) is also a polynomial in ar/(x) + b (a, b: real constants), we 
impose two condition (we assume G [xi,x 2 ]) 

rj(0) = and r]{x) : monotone increasing function, (2-19) 

which are not essential for (quasi-)exact solvability but important for expressing various 
formulas in a unified way. We impose another condition, to be called the symmetric shift- 
multiplication property: 

r]{x - i(3)r)(x + i/3) = (r](x) - rj(-i(3)) (r)(x) - r,(if3)) , (2.20) 

1 For the real shifts case, such r/(x) satisfying (|2.18[) and (|2.19[) can be classified into five types (|A.9p - (|A.13|) 
[5] and they also satisfy the condition (|2.20[) . 
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together with rj(x) ^ rj{x — i/3) ^ rj(x + i/3) ^ r](x). 

The two conditions (I2.18P and (I2.20p imply that any symmetric polynomial in rj(x — i/3) 
and r)(x + i/3) is expressed as a polynomial in r)(x). Especially we have (n > — 1) 

def rj(x — ij3) n+1 — r](x + i[3) n+1 ( a polynomial of" 
r](x — i/3) — 7](x + i/3) \degree n in r](x) 



9n{X) 



\n—k 



(2.21) 



fe=0 



The coefficient gn is real because g*(x) = g n (x). We set gh = except for < k < n. 
Since g n (x) satisfies the following three term recurrence relation 



9n+i{x) = (r)(x - i/3) + r)(x + i(3))g n (x) - r)(x - i/3)r](x + i(3)g n -i{x) (n > 0), 
we can write down g n explicitly Especially g^ for k — 0,1 are 



,(0) 



9 {l) 



[n + 1], 
±n(n + l)(2n+ l)ri 2 J 



for r\ 



(i) 



0, 



n[n+l] - (n + l)[n] j 2) m 



r_j for r} ^ 0. 



Here we have defined [n] as 



n 



def 



for = 0, 
for rj 1} > (< 



[ei — e ar 



>0)), 



for - 4 < rj 1J < («= ri ij = (e*2 - e^f ) 2 (0 < a < tt)). 

Note that and are expressed as 

r! 1} = [2] - 2, r (2 l = r]{-i(3) + r/(i/3). 

For n, m G Z, n > m — 1, we have 

^(ra + m + l)(n — m + l)(n 2 + 2n + m 2 ) 



for = 0, 



(i) 



(n + l)[n + 1] - m[m] - [|]- 2 [ !i± f ±i ][ ! ^f ±1 ] 
I0) 



r^i for 7^ 0. 



The following properties of [n] are useful: 



[a][a + c] - [6] [6 + c] = [a - 6] [a + fe + c] 

rn+m] rn— m+1 1 



Em 



[|] 



(n, m G Z, n > m — 1). 



(2.22) 

(2.23) 
(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 
(2.29) 
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2.3 potential functions 

The first goal is to construct a general form of the 'Hamiltonian' Ji such that a polynomial in 
r)(x) is mapped into another. It is achieved by the following form of the potential functions 

V±(x): 

V±(x) = v ±( x ) ^ ^2.30) 

(77(2; =f iff) — rj(x)) (r](x =F iff) — v( x =t 

V±{x)= VkMrfvixTiP) 1 , (2.31) 



fe,I>0 
k+l<L 



where L is a natural number roughly indicating the degree of rj(x) in V±(x) and v^i are real 
constants, with the constraint ^2 vf ;l ^ 0. It is important that the same Vj-j appears in 

k+l=L 

both V±(x). As we will see in the next subsection §2.4[ the 'Hamiltonian' 7i with the above 
V±(x) maps a degree n polynomial in r)(x) to a degree n + L — 2 polynomial (I2.36p . (12.391) . 

The essential part of the formula (12.301) is the denominators. They have the same form 
as the generic formula, derived by the present authors, for the coefficients of the three term 
recurrence relations of the orthogonal polynomials, (4.52) and (4.53) in [8j. The translation 
rules are the duality correspondence itself, (3.14)-(3.18) in [8]: 

S(n)^ V (x), -A n ^V + (x), -C n ^V-(x), 
a + (£(n)) — > rj(x — iff) — rj(x), a_(£(n)J — > rj(x + i/3) — rj(x). (2.32) 

Some of the parameters v k} i in (I2.3ip are redundant. From (12.181) and (12.201) . we have 

tj(x =F i/3) 2 = (2 + r { ^)7](x)r](x T ifi) - v( x ) 2 + r -l (v( x ) + v( x =F - Tl{-iP)rj{iP). (2.33) 

By using this repeatedly, a monomial r](x =F iff) 1 can be reduced to a polynomial of degree 
one in r](x =F iff) whose coefficients are polynomials in r)(x). Therefore it is sufficient to keep 
Vk,i with I = 0,1. The remaining 2L + 1 parameters Vkj (k + l < L, I = 0, 1) are independent, 
with one of which corresponds to the overall normalization of the Hamiltonian. In fact, if 
two sets of parameters {vk,i\ and {v' kl } (k + 1 < L, I = 0, 1) give the same V±(x), namely, 
Y,k=o( v k,o ~ v' kfi )r](x) k + Efc=oK,i ~ v 'k,i)v( x ) k v( x T W) = 0, then we obtain v kjl = v' k l . 
Therefore there is no more redundancy in v^i (k + 1 < L, I = 0, 1). Note that we have not 
yet imposed the boundary condition D(0) = (12.91) . The sinusoidal coordinate rj(x) itself 
may have extra parameters. 
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2.4 H on the polynomial space 

The action of U flZZEED on r]{x) n becomes with ( 12T30D and (12~3TD : 



^(x) n = e\V + (x) (r}(x - z/?) n - r](x) n ) + V-{x) (r)(x + i(3) n - r]{x) n ) 

V+{x) YTrZl V(x) r v(x - if3) n ~ l ~ r - V.{x) YT r Zl v(x) r v(x + i(3) n - l - r 



e 



rj(x — i(3) — rj(x + ij3) 



n-l 



Np Y V + (x)r](x - ifi) 71 - 1 - 7 ' - V_(x)r](x + ifi) 11 - 1 - 7 ' 



— — k v(x- i/3) l+n - 1 - r - r](x + i(3) l+n - 1 - r 

rj(x — i(3) — i](x + if3) 



E^T)(x) r E V k,l Vi x ) 



r=0 k,i>o 

k + l<L 



n-l 



k,l>0 r=0 
k+l<L 

= (a polynomial of degree n + L — 2 in r](x)) 

ri—l n+l—r—2 
k,l>0 r=0 j=0 

k + l<L 

n+L—2 m n—l 

= e e v(xr L - 2 - E E ^EfL-> 

m=0 j=max(m-L,0) M>0 r=0 

fc+Z — L — m+j 

n+L—2 m 

= E ^) n+i - 2 ~ m E <w.»- ( 2 - 34 ) 

m=0 j=max(m-L,0) 

Here e m Jin (the L-dependence is implicit) is defined by 

L—m+j n—l 



£ ^^-mEA-^' ( 2 - 35 ) 



i=0 r=0 

Therefore the matrix elements of 7i in the basis {i]{x) n } n= o t i„.. is given by 

n+L-2 n+L-2-m 

n v (xr= E v(x) m n\ n , n\ n = E e n+L - 2 _ mj ,„. (2.36) 

"1=0 j=max(n- 2— m,0) 

The coefficients e mj0i n and e m ^ n become, by using (I2.29P and (12 .2 7ft : 

jnj L-m 

e m ,o,n = Syfr E v L-m-i,i[ n+2 2 ~ L (2.37) 
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L— m+1 
Z=0 



(2) 
-1 



^n(n + 2/ - 2) ((n + / - l) 2 + I 2 - 2l)r 
+ + _ [I]-2[ 



for 



(i) 



0, 



ll-2rn+2f-2i rnl 



■2 J (2) 
' -1 



for 7^ 0. 



(2.38) 



So far the conditions = for / > 2 are not used. 
We have established 

iiV n C V n+L _ 2 , (2.39) 

where V n is the polynomial space defined in (12.171) . For L = 2, V n is 7Y- invariant. Therefore 
this case is exactly solvable; all the eigenvalues and eigenfunctions of 7i can be obtained 
explicitly and the eigenfunction is a polynomial of degree n in r](x) for each n. On the 
other hand, L > 3 cases are not exactly solvable but some cases can be made quasi-exactly 
solvable by certain modification to be discussed presently. For L = 0, 1 cases, the matrix 
o<m,n<K with finite K is not diagonalizable except for K = 0, 1. 
In the following we will set v^i = for I > 2, see §2.31 For the real shifts case, the 
condition D(0) = (12.91) is satisfied by choosing u ,o as v o,o — ^ v o,iV(~^)- 



3 Exactly Solvable 7i 

The L = 2 case is exactly solvable. Since the Hamiltonian of the polynomial space 7i is an 
upper triangular matrix f!2.36j) . its eigenvalues and eigenvectors are easily obtained explicitly, 
see Appendix O The eigenvalue £{n) is 



r n 1 



£(n)=m, 



e o,o . 



n-1 



[11 V^ 2 '° L - ; 



(3.1) 



and the corresponding eigenpolynomial P n (r](x)) is expressed as a determinant of the fol- 
lowing order n + 1 matrix, 

I w(x) r](x) 2 ■ ■ ■ r)(x) n 



PnHx)) 



oc 



£(0) - £{n) HI, 

£{l)-£{n) H\ 2 



-"0,2 



H v 
H v 



(3.2) 



£(n-l)-£(n) H^ 1>n 

For a choice of the sinusoidal coordinate among the possible forms flA.ip - flA.13P and the 
values of the five parameters, t> 0i o, ^i.Oj ^o.i' ^1,1 an d t '2,o, these two formulas (13.11) and (13. 2p . 
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although clumsy, give the complete solutions of the 'Schrodinger equation' TtP n (r](x)) = 
£(ri)P n (r)(x)) at the algebraic level. For the solutions of a full quantum mechanical problem, 
however, one needs the square-integrable groundstate wavefunction <po(x) (12.14j) . which is 
essential for the existence of the Hamiltonian TC and the verification of its hermiticity. These 
conditions would usually restrict the ranges of the parameters i> 0) o, • • • , ^2,0- 

For specific problems, however, there are more powerful and systematic solution methods 
based on the shape invariance [H [6J, [TJ, [8J, [9] and the closure relation [5j [8j [9] . These two are 
independent and sufficient conditions for exact solvability which are applicable to not only 
ordinary QM but also discrete QM. In our previous works [61 [TJ El [8j [9] these conditions were 
verified for each specific problem. Here we will provide proofs based on the generic form of 
the exactly solvable (L = 2) 'Hamiltonian' TC, (I2.16p . (I2.30p . (12.311) . These proofs apply to 
all the exactly solvable discrete QM. In the rest of this section we assume the existence of 
the ground state wavefunction <fio(x) (12.141) . 

3.1 closure relation 

The closure relation is a commutator relation between the Hamiltonian TC and the sinusoidal 
coordinate r](x) 

[TC, [H,rj\) = r)Ra(H) + [H,rj\ R\{TC) + R^(TC). (3.3) 

Here R\(z) are polynomials with real coefficients r!f\ 

R^z) = r?z + r{°\ R (z) = rfz 2 + z + rj 0) , R^(z) = r^z 2 + r^z + r^. (3.4) 

Reflecting the fact that the Hamiltonian TC has shift operators e ±/3p , whereas rj(x) has none, 
the function Ro{7i) and R-iiTi) are quadratic in TL and R\{7i) is linear in Ti. By similarity 
transforming (13.31) in terms of the ground state wavefunction 0o> h is rewritten as 

[H,[H,r ] ]] = v R (H) + [H,r ] ]R 1 (H) + R- 1 (H). (3.5) 

The closure relation (13. 3p allows us to obtain the exact Heisenberg operator solution for rj(x), 
and the annihilation and creation operators are extracted from this exact Heisenberg 
operator solution [5]: 

e im 7](x)e- im = a^e ia+{n)t + a { - ] e ia - Wt - R^H) RoiH)- 1 , (3.6) 
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a±(H) = \(R X (K) ± y/Ri(H) 2 + 4R (H) ), (3.7) 

Ri(H) = a + (H) + a^H), R (H) = -a + {H)ot-{H), (3.8) 

a (±) drf ±^ n ^( x ^ _ [ v (x) + R^i^RoiHy^a^H)^ (a+(n) - a^H))' 1 (3.9) 

= ±(a+(H) - a-{H))~ l ([H^x)] +a ± {H){i 1 {x) + R- 1 {H)R {Hy 1 )}. (3.10) 



The energy spectrum is determined by the over-determined recursion relations £{n + 1) = 
S(n) + a + (S(n)) and £{n — 1) = £ (n) + a_(£(n)) with £(0) = 0, and the excited state 
wavefunctions {(b n (x)} are obtained by successive action of the creation operator on the 
groundstate wavefunction 4>q(x). The closure relation (13.51) (or (I3.3P ) is equivalent to the 
following set of five equations: 

rj(x — 2i/3) — 2rj(x — i/3) + 77(a) = r^r)(x) + + (n(x — i/3) — 77(a)) , (3-11) 
i](x + 2i/3) — 2r](x + i/3) + r](x) = 77(a) + + (r](x + i/3) — 77(a)) , (3.12) 
(rj(x - i/3) - 77(a)) (V+(x - i/3) + V-(x - i/3) - V+(x) - V-(x)) 
= -(r (2) r](x) + r {2 \) (V + (x - i(3) + V.{x - i/3) + V+(x) + V-(x)) 

- (r](x - %I3) - 77(a)) (V+(x - i(3) + V-(x - i/3)) 

+ e- 1 (r { 1) r ] (x) + r a \ + rf } (r/(a - i/3) -77(a))), (3.13) 
(r](x + i/3) - r](x)) (V+(x + i/3) + V-(x + i/3) - V+(x) - VL(x)) 
= -(r®r)(x) + r { $) (V+(x + i/3) + V^{x + i/3) + V+(x) + V_{x)) 

- (r](x + i/3) - 77(a)) (V+(x + i(3) + V-(x + i/3)) 

+ e- 1 (r ( Q ) r ] (x) + r { }\ + rj 0) (r/(a + i/3) -77(a))), (3.14) 
2(?7(x) - ?7(x - i/3)) F+(a)V_(a - i/3) + 2 (77(a) - r/(a + i/3))y_(a)V+(a + i/3) 
= (rf ^(x) + rf 2 )) (V+(x)V_(x - i/3) + VL(x)V + (x + i/3) + (V + (z) + V^x)) 2 ) 
+ rj 1} (?7(x - i/3) - 77(x))V+(x)K-(x - i/3) + rf ] [t]{x + i/3) - 77(a)) V_(a)V+(a + if3) 

- e- l {r$\(x) + r3) (V + (a) + K-(x)) + e~ 2 : (r£ o: 1 77(a) + rl°J). (3.15) 

Obviously (13.111) and (13.121) are equivalent and so are (I3.13P and (13 . 14j) . under the condi- 
tion (I3.16p . By substituting our choice of V±(x) (I2.30j) (I2.31j) for L = 2, it is straightforward 
to verify the other three equations (I3.13p - (l3.15p . The coefficients appearing in (13.41) are 
expressed by the parameters vx,q, fo,i, and t>2,o together with the two parameters r± 
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and r_{ which have already appeared in the definition of r)(x) (12. 18ft (see also ( 12.261) ): 

rj 1) = 2rf ) , (3.16) 
e~ 2 r { Q ] = -v 2 ,qVi j1 , (3.17) 
e-V°l = -v 2 , v 0>1 . (3.18) 

Note that vo,o does not appear. It implies that for the imaginary shifts case the commutation 
relation between the annihilation and creation operators does not depend on vo,o- With these 
formulas, the explicit forms of a±(7i) (13. 7p can be expressed in terms of \ t^o and v^x. 
It is straightforward to verify the eigenvalue formula (13. ip . This concludes the unified proof 
of the closure relation for all the discrete QM. 

3.2 dual closure relation 

The dual closure relation has the same forms as the closure relation (13.31) and (13. 5p with the 
roles of Hamiltonian 7i (71) and the sinusoidal coordinate rj(x) exchanged: 

[77, [77, H]]=H < ial (r/) + [77, H\ < al (r/) + R d ™\ V ), (3.19) 
[77, [77, H\\=H < ial (77) + [77, H\ R d ^(rj) + R d ™\ V ), (3.20) 

where R dnal (z) are as yet unknown polynomials. We will show below that the dual closure 
relation is the characteristic feature shared by all the 'Hamiltonians' 7i which map a poly- 
nomial in 77(2;) into another. Therefore its dynamical contents are not so constraining as the 
closure relation, except for the real shifts (the discrete variable) exactly solvable (L = 2) 
case, where the closure relation and the dual closure relations are on the same footing as 
shown in [§]. By substituting the 'Hamiltonian' 7i (12.161) without any further specification 
of V± into the above (I3.20p . we find it is equivalent to the following set of three equations: 

(r)(x) - r](x - 7/3)) 2 = R dnal (r,(x - i(3)) + (r,(x) - 77(0; - ip))R A ™ l {r,(x - i/3)), (3.21) 
(r,(x) - v (x + if3)) 2 = R dnal (r,(x + i(5)) + (rj(x) - ^x + i(3))R dnal (r,(x + if3)), (3.22) 
= -e(V+(x) + y„(x))^ ual (77(a;)) + R d _ u f (r]{x)) . (3.23) 

These imply 



,(2) 



,(1) 

1 > 



-U°) 



^2,0 + Ul,l> 



-1 (1) 1 
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#o ual 0?(X>) =-(-n(x-ip)-Ti(x))(Ti(x + il3)-ri(x)), (3.25) 
R*f(r,(x)) = e{V + (x) + V.(x))R^\ V (x)). (3.26) 

By using the denning properties of the sinusoidal coordinate ( I2.18I) - (I2.20I) . we actually find 
that Rf nal (z) is a degree 1 polynomial in z and i?Q Ual (z) is a quadratic polynomial: 

R<t™\z) =rPz + r { *l, (3.27) 
R^\z) = rjV + 2r m lZ - r){-ip)r){ip). (3.28) 

By using the explicit forms of V± (l2.3U|) - (j2.31j) (with an arbitrary L) we obtain 

L 

R*f( z ) = e(y , + ^2(v kfl + (3.29) 
fc=i 

Therefore all Rf nal (z) are polynomials and the dual closure relation is demonstrated in a 
unified fashion for an arbitrary L. Thus it does not characterise the exact nor the quasi- 
exact solvability. For the exactly solvable L = 2 case, by using ( 13. 17ft and ( I3.18p . R < zf 1 (z) 
can be written as 

B*f(z) = rf ] z 2 + r W lZ + ev 0jQ . (3.30) 

For the real shifts case, in order to satisfy D(0) = (12.91) . we have to take f 0) o = — T](— l)fo,i 
and this implies ev 0fi = r)(l)ri{-l)B(0). See (4.104)-(4.106) in [5]. 

3.3 Askey- Wilson algebra 

Here we will focus on the exactly solvable systems and will briefly comment on the rela- 
tionship between the closure plus the dual closure relations and the so-called Askey- Wilson 
algebra [TSl dSl EDI EI] • By simply expanding the double commutators in the closure (13.31) 
and the dual closure (I3.19P relations, we obtain two cubic relations generated by the two 
operators H and rj: 

H 2 r] - (2 + rP)HvH + r(H 2 - rf\Uq + rfH) - rj% = r {2 \H 2 + r^H + r ( °\, (3.31) 

rfU - (2 + rP)r]Hrj + Hrf - r^H + Hrf) + rj(-i(3)ri(i/3)H = r?V + r^q + ev 0fi . 

(3.32) 

From its structure, the closure relation is at most linear in r\ and at most quadratic in 7i. 
So the l.h.s. of (I3.3ip has terms containing one factor of rj and the r.h.s, none. It is simply 
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R-i(7i). Likewise, the l.h.s. of (13.321) has terms containing one factor of Ti and the r.h.s, 
none. It is simply R^ x l (Tl). In (13.321) . r){—i(3)r]{i(3) is just a real number, not an operator. 

These have the same form as the so-called Askey- Wilson algebra, which has many different 
expressions. The original one is due to Zhedanov [18] . Here we present a slightly more general 
version than the original one and is due to [T9], [20] . It is generated by three elements K~i, 
K 2 , K 3 : 

[K 1 ,K 2 }=K 3 , (3.33) 
[K 3 , K x ] = 2pK 1 K 2 K 1 + a 2 (K x K 2 + K 2 K X ) + a x K 2 + c 2 K 2 + dK x + g 2 , (3.34) 
[K 2 , K 3 \ = 2pK 2 K 1 K 2 + a x (K 2 K x + K X K 2 ) + a 2 K\ + c x K x + dK 2 + g x . (3.35) 

By expanding the commutators and eliminating K 3 , they are reduced to 

K\K 2 + 2(p - l)K x K 2 K x + K 2 K\ + a 2 {K x K 2 + K 2 K X ) + c 2 K 2 = -a x K\ - dK x - g 2 , (3.36) 
K\K X + 2(p - l)K 2 K x K 2 + K X K\ + a x {K 2 K x + K X K 2 ) + c x K x = -a 2 K 2 2 - dK 2 - g x . (3.37) 

Another version due to Terwilliger [21] is generated by two independent elements A and A x 
and it has only expanded forms: 

A 2 A X - p T AA x A + A X A 2 - •y(AA x + A x A) - p T A x = 7 M 2 + uA + r) T , (3.38) 
A x2 A-(3 T A x AA x +AA x2 - 1 x (A x A + AA x )-p x A = 1 A x2 + uA x +ri x . (3.39) 

Here is the list of correspondence of the generators and coefficients: 



ref . PH ED] 


ref. [21] 


this paper 


K x 


A 


n 


K 2 


A x 


V 


2(1 -p) 


Pt 


2 + rf } 


—a 2 


7 


rf ] = e(v 2fi + v XjX ) 


—a x 


7 X 


r {2 l = V (-i(3) + V (i(3) 

(0) 2 


-c 2 


Pt 


-c x 


Pt 


-r](-i(3)r](i/3) 


-d 


uJ 


r_ x = e(v Xfi + v QjX ) 

(0) 2 




Vt 


-9i 


Vt 





(3.40) 



In rrj] the Casimir operator Q commuting with all the generators of the algebra, [K x , Q] = 
[K 2 ,Q] = [K 3 ,Q]=0is given: 

Q = K X K 2 K X K 2 + K 2 K X K 2 K X - (1 - p)(K x K 2 2 K x + K 2 K 2 X K 2 ) 
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+ (2 - p)(a 1 K 1 K 2 K 1 + a 2 K 2 K 1 K 2 ) + (1 - p)(c x K\ + c 2 K\) (3.41) 
+ (d- a x a 2 ){K x K 2 + K 2 K X ) + ((2 - p) gi - a 2 c 1 )K 1 + ((2 - p)g 2 - a x c 2 )K 2 . 



With the above substitution (13.401) . K\ — > H, K 2 — > 77, etc, the Casimir operator turns out 
to be a constant [22] : 

Q = e 2 (wi,ifo,o - v h0 v ,i - r^lv 2 ,oV 0jl ). (3.42) 

Although this fact might appear striking from the pure algebra point of view (I3.33p - (l3.35p . 
it is rather trivial in quantum mechanics. In one- dimensional quantum mechanics, there is 
no dynamical operator which commutes with the Hamiltonian. Therefore, if Q commutes 
with Tt, it must be a constant. 

Now here are some comments on the dissimilarity. The first obvious difference is the 
structure. While the Askey- Wilson algebra f!3.33j) (I3.35j) or (I3.38I) - (I3.39I) has no inherent 
structure, the closure relation (I3.3P has the right structure to lead to the Heisenberg operator 
solution for r](x), whose positive and negative energy parts are the annihilation-creation 
operators [8], [9]. It is the Hamiltonian and the annihilation-creation operators that form 
the dynamical symmetry algebra of the system [HIS], not the closure or dual-closure relations, 
nor the Askey- Wilson algebra relations. The g-oscillator algebra of [23] is the typical example 
of the dynamical symmetry algebra thus obtained. 

The next is the difference in character of the Askey- Wilson algebra itself for the two 
cases; the pure imaginary shifts and the real shifts cases. The main scene of application of 
the Askey- Wilson algebra is the theory of the orthogonal polynomials of a discrete variable. 
The (g-)Racah polynomials are the typical example of this group [131 E]- In our language, it 
is the theory of the eigenpolynomials of 7i in discrete quantum mechanics with real shifts. 
One outstanding feature of these polynomials is the duality [SI EI]. For the eigenpolynomials 
of U 

HP n (v(x)) = £(n)P n ( V (x)), n = 0, 1, ... , (3.43) 
there exist the dual polynomials Q x (£(n)), satisfying the relation 

P n {7i(x)) = Q x {e{n)), x = 0,l,..., n = 0,l,.... (3.44) 

This duality x <-> n, rj ~ r]{x) <-> £ (n) ~ Ti is reflected in the symmetry between the pair of 
operators (called the Leonard pair [24J) K\ and K 2 or A and A x in the Askey- Wilson algebra. 
The Askey- Wilson algebra or the closure and dual closure relations are quite instrumental 
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in clarifying various properties of the pair of orthogonal polynomials of a discrete variable 



Now let us consider the discrete quantum mechanics with the pure imaginary shifts. 
In this case, the sinusoidal coordinate rj(x) takes the continuous value (spectrum) for the 
continuous range of x G (xi,^) (12.21) . which is markedly different from the spectrum of H, 
postulated to take the semi-infinite discrete values (12.61) . The eigenpolynomials TCP n (r](x)) = 
£(n)P n (r)(x)) depend on the continuous parameter x and they have no dual polynomials. 
The Askey- Wilson and the Wilson polynomials are the typical examples [131 IS] ■ As shown 
in previous work [51 [9] and in §3.21 the essential information on exact solvability is contained 
only in the closure relation (13. 3p . There is no evidence that the dual closure relation plays 
a comparable role to the closure relation. Therefore we may conclude that the apparent 
symmetry between Ti and rj, or K\ and K 2 or A and A x in the Askey- Wilson algebra is 
quite misleading for the pure imaginary shifts case. In other words, a part of the Askey- 
Wilson algebra is irrelevant to the orthogonal polynomials of a continuous variable. 

3.4 shape invar iance 

Let us briefly review the condition and the outcome of the shape invariance [1] in our lan- 
guage. In many cases the Hamiltonian contains some parameter (s), A = (Ai, A2, • • •)• Here 
we write parameter dependence explicitly, 7i(A), ^4(A), £ {n ; A), (f> n (x ; A), etc, since it is the 
central issue. The shape invariance condition is [H [7J [8j [9] 



where k is a real positive parameter and A' is uniquely determined by A. Let us write the 
mapping as a function, A' = si(A). In concrete examples, if we take A appropriately, A' has 
a simple additive form A' = A + 8. The energy spectrum and the excited state wavefunction 
are determined by the data of the groundstate wavefunction (f) (x A) and the energy of the 
first excited state £{1 ; A) as follows: 



mm- 



kA(\'YA(\')+£(1;\) 



(3.45) 



n-1 




(3.46) 



s=0 



0„(x ; A) oc ^(A [0] )U(A [11 )U(A [2] ) t • • ■ ^(A [n - 1] )Vo(^ ; A 



[n] 



)• 



(3.47) 



Here A N is A 



A, A N =si(A [n ~ 1] ) (n = 1,2,...). 
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3.4.1 pure imaginary shifts case 

Here is a unified proof of the shape invariance for the discrete quantum mechanics with pure 
imaginary shifts. The shape invariance condition (I3.45P is decomposed to the following set 
of two equations: 

V(x - q ; X)V*(x -q ] \) = K 2 V(x ; X')V*(x - ij ; A'), (3.48) 
V{x + q-,X) + V*(x — q ■ A) = k{V(x ■ A') + V*{x ; A')) - S(l ; A). (3.49) 

We assume that r](x) satisfies the relation 

rj(x) = [~] (r](x - q) + r](x + q) - r](—q) - r](q)) . (3.50) 

Moreover we assume that r)(x; A') = r](x;X), for example it is satisfied if rj(x) is A- 
independent. Both are easily verified in each of the explicit examples listed in the Ap- 
pendix El flA~T]) - flA~8|) . When the forms of the potential functions (12301) and (HOB (with 
L = 2) are substituted, the shape invariance conditions (I3.48I) - (I3.49I) are satisfied. If we take 
{vk,o (k = 0, 1, 2), Vk,i {k = 0, 1)} as A, then A' and £{1 ; A) are 

KV 2fl = -vi,u ( 3 - 51 ) 
«<i = «a,o + [2Ki, (3.52) 

' " [i]Ko-«b,i)+*--}f^T«a.o + «i,i). ( 3 - 53 ) 



KV 



1.0 



L 2 J 

ii2 r 1 1 r 3 1 



<i = [§Ro + [§H,i + rS ^2,0 + ^ «i,i) . ( 3 - 54 ) 

/ , (2)/[|][|] [l] 2 \ . rH2/[l] 4 (2)2 / \ / - \ \ 

Kv o,o = ^o,o + rl{ [-rrr «o,i - Tnr v i,o J + [5] (^TTu r -i ~ V{-H)v(W J ^2,0 



L 2 J L 2 J L2 

■11 / L4H4J J2)2 . r 3 



L2H [IJ3 



^S 2 + [|]»7(-i7)»?(i7))vi,i s ( 3 - 55 ) 



£(l;A) = vi,i. (3.56) 

Note that the above formula £ (1 ; A) is consistent with the general formula (13.11) . It is 
elementary to verify that the quadratic recursion formula generated by (I3.5ip and (13.521) 
coupled with the shape invariance energy formulas (13.461) and (I3.56P reproduces the energy 
eigenvalue formula (13.11) . However, the other formulas for the parameter shifts (I3.53I) - (I3.55I) 
seem too complicated to be practical. As shown in 0, [3, [5], [9] , the parameter shifts are much 
simpler for the known examples. 
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remark In ordinary quantum mechanics there is a method for constructing a family of 
isospectral Hamiltonians, known as Crum's theorem [2]. Recently we have obtained its 
discrete quantum mechanics version, see [25]. If (fii(x) take a form </>i(x) = (fio{x) (const + 
const ■ i](x)), which occurs indeed in the setting of this paper, the potential function of the 
first associated Hamiltonian is given by 



V 1 Hx + ij) = V{x) — — — 

7]{X) — Tj\X + Z7J 



(3.57) 



Therefore, if shape invariance holds, V(x) satisfies 

in which the sinusoidal coordinate may depend on A. 
3.4.2 real shifts case 

The shape invariance (13.451) is equivalent to the following set of two equations: 

B{x + 1; X)D{x + 1 ; A) = k 2 B(x ; \')D{x + 1 ; A'), (3.59) 

B(x ; A) + D(x + 1 ; A) = k(B(x ; A') + D(x ; A')) + £(l ; A). (3.60) 

For the classified five types of r)(x), (i)'-(v)' in (1A.9|) - (1A.13I) . the shape invariance holds. 



The boundary condition D(0) = (12.91) forces to choose v 0t0 as v 0t0 = —v 0) ir](—l). Thus we 
take the parameters {vk$ {k = 1,2), (A; = 0, 1)} (and d for (ii)' and (v)') as A, then A' 
and £{1 ; A) are 

KV 2,o = — fi,i, (3.61) 
Ki=^ + [2Ki, (3.62) 



KV 



/"[l](^i,o - ^o,i) + V[\\v{ l )v2fi + vr^v^i, (3.63) 



1,0 — PL^JV^O U 0,1J ~ PL2 



«<i = M([fHo + [§Ki) + M[|]^(l)^2,o + (^rLl + /i[§] (77(1) - t/C-I)))^!,!, (3.64) 
f (1 ; A) = (3.65) 



d + 1 for (ii)' 



|<ig for (v)', ^ ^ 

in which and v are constants 

1 for (i)'-(ii)' ( x for (i)'-(iv)' 

fi = { g-a for (ill)' 1/ = <j 1 + for , (3.67) 



for (iv) - (v)', Kl + d 
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The quadratic recursion formula (13.611) and f!3.62j) are exactly the same as those of the pure 
imaginary shifts case (I3.5ip and (13.521) . Therefore the shape invariance energy formulas 
(I3.46P and (13.651) produce the same energy spectra (13.11) . For the finite dimensional case, 
the natural number TV satisfying B(N ; A) = (12.91) is also counted as a varying parameter. 
Then the shape invariance including the conditions 

N' = N-1, B(N';X') = (3.68) 

is satisfied. As shown in [8] , the parameter shifts are much simpler for the known examples. 

4 Quasi-Exactly Solvable TC 

Quasi-exact solvability (QES) means that only a finite part of the spectrum and the cor- 
responding eigenfunctions can be obtained exactly |15j . Usually such a theory contains a 
finite dimensional vector space [26J consisting of polynomials of a certain degree which forms 
an invariant subspace of the 'Hamiltonian' Tl, or more precisely its modification TC . There 
are many ways to accomplish QES. The method of this paper can be considered as a simple 
generalisation of the one in [27]. That is, to add non-solvable higher order term(s) together 
with compensation term(s) to an exactly solvable theory. As is clear from the construction, 
the sinusoidal coordinate plays an essential role. 

For a given positive integer M, let us try to find a QES 'Hamiltonian' TC, or more precisely 
its modification TC ; , having an invariant polynomial space Vm'- 

H'V m C V m - (4.1) 

For L > 3, (127341) is 

L— 3 m 

Hi](x) n = ^2 v( x ) n+L ~ 2 ~ m ^2 em 'i> n + ( a P°ly nomia l of degree n in rj(x)) . (4.2) 

m=0 j=0 

So let us define TC' by adding compensation terms to TC as 

L— 3 m 

n><^H-J2e m (M)r)(x) L - 2 - m , e m (M) d ^ ^ e md , M . (4.3) 

m=0 j=0 

Then we have TC'r](x) M G Vm- For 1 < m' < L — 3, we have 

L—m'—3 m 

H' V (x) M ~ m '= Yl v{x) M+L - m '- 2 - m (^e m , hM . ml -e m {M)) 

m=0 j=0 
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+ (a polynomial of degree M in T}(x)). (4-4) 
If we could choose Vk t i to satisfy all these conditions 

m 

e m,j,M-m' - e m {M) = (1 < m' < L - 3, < m < L - m' - 3), (4.5) 

j=0 

then we would obtain TL'Vm Q Vm- 

4.1 QES with L = 3 

For the L = 3 case, W is defined by adding one compensation term of degree one 

H> d ^H- e (M) v (x), eo(M) = e , ,M, (4.6) 

and we have achieved the quasi-exact solvability WVm Q Vm- The number of exactly 
determined eigenstates is M + 1 = dim(Vjvf)- in this case there is no extra conditions for 
Vk,i- The explicit form of eo(M) is 

eo(M)= £ l|l([^Ko+[^Ki). (4.7) 

This QES theory has two more parameters i^o and i>2,i on top of those in the original 
exactly-solvable theory (L = 2). Most known examples of QES belong to this category but 
those in ordinary quantum mechanics have only one extra parameter. 

4.2 QES with L = 4 

This type of QES theory is new. For L = 4 case, W is defined by adding a linear and a 
quadratic in rj(x) compensation terms to the Hamiltonian 7i: 

H' d ^H- e (M) v (x) 2 - ei (M)r)(x), e (M) = e , ,M, e 1 {M) = e 1AM + e ljl)M) (4.8) 

and H'i](x) M G V M . By using fl2T28|) we have 

H'r)(x) M ~ l = rj(x) M+1 (e ,o,M-i — e (M)) + (a polynomial of degree M in r](x)^j 

= -Sf](x) M+1 ([M - l}v 4 ,o + [M]v 3A ) + (a polynomial of degree M in r](x)). (4.9) 

In order to eliminate the r](x) M+l term, we choose ifyi as 

\M — 11 . . 

w 3 , 1 = - L -p^f4,o. (4-10) 
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We have achieved the quasi-exact solvability H'Vm Q Vm- The explicit forms of eo(M) and 
ei(M) are 

e (M) = -e \ r 2J[ 2 ^ 40 , (4.11) 
oi ; [|][M + 3] ' V ; 

eiW = e^([^Ko + [^ 1 M 

A/(A/-l)(Af 2 +5Af+8) r (1) _ n 

(2) I Af+3 1 ~~ ' 

Er-1«4,0 X { 2[«][ifcl] „.„, , nrl ,[2M + 5]\ . (1) , „ (4-12) 



i S^ l [4]-2[3] + 2[l 1 pp) forr'Vo. 



The theory has three more free parameters on top of those of the original exactly solvable 
theory (L = 2). 

4.3 non-QES for L > 5 

The higher L becomes, the number of conditions to be satisfied (14. 5p increases more rapidly 
than the number of additional parameters. We will show that L > 5 case cannot be made 
QES. The condition (14. 5 j) with m = gives eo,o,Af = e,a,<o,M-m' (1 < m' < L — 3), and by 
using (12371) and (12T28|) we obtain 

[ M - ^Ko + [M - ^H-i,i = (l<m'<L-3). (4.13) 

For L > 5 case, these equations do not have non-trivial solutions. For m' = 1,2 we obtain 

[M - 1] [M] \ ( v Lfi \ _ fO 



[M - \\) \0J- (4 ' 14) 

The determinant of this matrix is [|] which does not vanish. Thus we obtain vt,o = = 
0. Namely there is no Vk,i (k + l = L) term. Therefore L > 5 case cannot be made QES. 

5 (Quasi-) Exactly Solvable Hamiltonian 

If there exists a groundstate wavefunction <f>o(x) which satisfies (12.141) (and ||0o|| < oo, the 
hermiticity of TC), we can return to the Hamiltonian TC from the 'Hamiltonian' TC by the 
inverse similarity transformation (12.151) . In the same way the QES Hamiltonian TC is ob- 
tained from TC by the inverse similarity transformation in terms of the pseudo- groundstate 
wavefunction <f>o(x) satisfying A(f)o = (12.141) . 

H'^^oW'o^)- 1 . (5.1) 
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It should be noted that 4>o{x) is neither the groundstate nor an eigenstate of the total 
Hamiltonian Ti' . Thus it is called the pseudo-groundstate wavefunction. For the L = 3,4 
cases we have 



L = 3 : H = H- e (M)r)(x), 

L = 4 : H' = H-e (M)r](x) 2 - ei (M)r](x). 



(5.2) 
(5.3) 



Let us note that the Hamiltonian W does not factorise and the semi positive-definite spec- 
trum is lost due to the compensation terms. For the pure imaginary shifts case, the existence 
of (pseudo-)groundstate wavefunction 4>o(x) strongly depends on the concrete form of V(x) 
and its parameter range. There is no general formula to write down 4>q(x) in terms of V(x). 
On the other hand, for the real shifts case, the (pseudo-)groundstate wavefunction 4>o(x) is 
uniquely given by [8] 



The positivity of B(x) and D(x) (\2.9\\ restricts their parameter range. For the infinite case 
x G [0,oo), the square- summability ||0o|| < oo restricts the asymptotic forms of B(x) and 



6 Summary and the Recipe 

Based on the sinusoidal coordinate rj(x) we have systematically explored a unified theory of 
one-dimensional exactly and quasi-exactly solvable 'discrete' quantum mechanical models. 
The Hamiltonians of discrete quantum mechanics have shift operators as exponentiated forms 
of the momentum operator p = —id x , e ±l3p = e Tll39x . This method applies to both the pure 
imaginary shifts (j3 = 7 6 K^o) an d the real shifts cases (/3 = i = yj— 1), which have a 
continuous and a discrete dynamical variable x, respectively. The main input is the special 
form of the potential functions V±(x) f)2.30p and (12.311) . with which the 'Hamiltonian' 7i 
(I2.16P maps a polynomial in r/(x) into another. We obtain exactly solvable models (degree 
n — > n) and quasi-exactly solvable models (degree n-^n+l,n-|-2)by adding compensation 
terms, which are linear and quadratic in f](x), respectively. The QES Hamiltonians based 
on the mapping (degree n — > n + 2) are new. 

The corresponding result in the ordinary QM can be found in the Appendix of [5J. This 



x-l 



B(y) 




(5.4) 
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early work, however, does not cover the quasi-exactly solvable cases. In this connection, see 
the recent developments [28] . 

The present paper is for the presentation of the basic formalism. Application and concrete 
examples will be explored in a subsequent publication [H]. The explicit forms of various 
sinusoidal coordinates are listed in Appendix lAl 

The forms of the (pseudo-)groundstate wavefunctions (j)o(x) for the pure imaginary shifts 
(the continuous variable) case depend on the choice of the sinusoidal coordinates. They are 
'gamma functions' having various shift properties; the (Euler) gamma function for (i)-(ii) 
(IA.lj) - (lA.2p . the (/-gamma function (or the g-Pochhammer symbol) for (iii)-(iv) (1A.3|) - (1A.4|) . 



the double gamma function (or the quantum dilogarithm function) for (v)-(viii) flA.5j) - flA.8j) 



In a subsequent publication we will present explicit examples of new Hamiltonians based on 
flA.5j) - flA.8j) . Their eigenfunctions contain orthogonal polynomials and the double gamma 
functions as the orthogonality measure functions. Eigenpolynomials P n (r](x)) for exactly 
solvable QM belong to the Askey-scheme of hypergeometric orthogonal polynomials. Various 
examples of exactly solvable QM were investigated for (i)-(iv) (1A.lj) - (1A.4l) |9J and (i) — (v) 
flA.9j) - flA.13j) [8J, and quasi-exactly solvable QM were partially examined in [T6| IT?]. 

The simple recipe to construct an exactly or quasi-exactly solvable Hamiltonian is as 
follows: 



'1) Choose the sinusoidal coordinate among (1A.lj) - (lA.8l) if the variable is continuous 



among flA.9l) - flA.13j) for the discrete variable. 



(2) Choose L = 2 for exact solvability and L = 3,4 for quasi-exact solvability and write 
down the 'Hamiltonian' 7i in the polynomial space (12.301) and (I2.3ip with the free 
parameters v^i, k + 1 < L, I — 0, 1. For the quasi-exactly solvable case, add the proper 
compensation terms (14. 6 1) - (14. 7p for the L = 3 case and (I4.8I) - (I4.12I) for the L = 4 case. 

(3) Determine the (pseudo-)groundstate 0o as a zero mode of A, (I2.14p . which can be found 
among the various gamma functions listed as above or (15.41) for the discrete variable 
case. 

(4) Restrict the parameter ranges so that the square-integrability of (fio and the hermiticity 
is satisfied for the continuous variable case. For the discrete variable case the positivity 
B(x) > and D(x) > and the boundary condition(s) D(0) = 0, (B(N) = 0) (T23D 
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are the conditions to restrict the parameters. For the infinite dimensional case the 
square summability J2T=o fioiz) < 00 mus t be satisfied, too. 

(5) Apply the inverse similarity transformation (15. ip in terms of 0o 011 the 'Hamiltonian' 
in the polynomial space to get the Hamiltonian 7i or TC'. 

Acknowledgements 

We thank Paul Terwilliger for fruitful discussion. This work is supported in part by Grants- 
in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and 
Technology, No.18340061 and No.19540179. 

A List of Sinusoidal Coordinates 

Here we list the explicit forms of the sinusoidal coordinates, based on which various concrete 
examples of exactly and quasi-exactly solvable theories are constructed. 

These are eight sinusoidal coordinates for the pure imaginary shifts case (continuous x): 



(i) : T)(x) = x, —00 <x < 00, 7 = 1, (A.l) 

(ii) : i](x) = x 2 , 0<x<oo, 7 = 1, (A. 2) 

(iii) : rj(x) = 1 — cosx, <x < 7r, 7 G K^o, (A-3) 

(iv) : r](x) = smx, -§<cc<f, 7 e R^ , (A.4) 

(v) : r](x) = l-e~ x , -oo<x<oo, 7 G E^ 0) (A. 5) 

(vi) : rj(x) = e x - 1, -00 < x < 00, 7 G R^ , (A. 6) 

(vii) : T)(x) = coshx — 1, 0<x<oo, 7G M^o ; (A. 7) 

(viii) : r)(x) = smhx, —00 <x < 00, 7 G BL^cb (^--8) 

and five sinusoidal coordinates for the real shifts case (integer x): 

(i) ' : r](x) = x, (A.9) 

, , r 1 for d > — 1, 

(ii) ' : rj(x) = e'x(x + d), e' = \ (A.10) 
V ; /V ; V ' I -1 for d < -N, V ' 

(iii) ': r]{x) = l-q x , (A.ll) 

(iv) ': rj(x) = q~ x - 1, (A.12) 
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v 



r 1 for d < q 1 , 
.': V (x)=e'(q~*-l)(l-dq*), e> = { \^ (A.13) 

l> —1 tor a > (7 . 



1 for d > q 



where < q < 1. As shown in detail in §4C of [8], the above five sinusoidal coordinates for the 
real shifts (IA.9p - flA.13p exhaust all the solutions of fl2.18p - fl2.20p up to a multiplicative factor. 



On the other hand, those for the pure imaginary shifts (i)-(viii) (1A. lh — (1A.8I) are merely 
typical examples satisfying all the postulates for the sinusoidal coordinate fl2.18p - fl2.20p and 
the extra one used for the shape invariance (13.501) . It is easy to see that i](x) = x + sinh(27rx), 
— oo < x < oo, 7=1 is a good sinusoidal coordinate for the imaginary shifts but it fails to 
fulfill the extra condition fl3.50p . 

B Hermiticity of the Hamiltonian 

In this Appendix we recapitulate the proof of the hermiticity of the Hamiltonian 7i ( 12. 71) 
for the pure imaginary shifts (continuous variable) case [T6l 19]. For the real shifts (discrete 
variable) case the Hamiltonian 7i (12.71) is a hermitian matrix (real symmetric matrix) and 
there is no problem for the hermiticity. Thus we consider only the continuous variables 
case, in which the wavefunctions and the potential functions are analytic function of x as 
explained in §2J The *-operation on analytic functions is also defined there. 

By using the formula (AB)^ = B*A\ we obtain Ti) = Ti but this is formal hermiticity. 
In order to demonstrate the true hermiticity of H, we have to show (g,T~Cf) = (Hg,f) 
with respect to the inner product (12.41) . Since the eigenfunctions considered in this paper 
have the form 4>o(x)P n (r](x)) , it is sufficient to check for f(x) = <p (x)P(r](x)) and g(x) = 
4>o{ x )Q(v{ x )) j where P{rj) and Q(rj) are polynomials in 77. Since H is real H* = H, namely 
TC maps a 'real' function to a 'real' function (f*(x) = f(x) =3- (Hf)* = (Hf)), we can take 
4>o(x), P(rj(x)) and Q{j]{x)\ to be 'real' functions of x and we do so in the following. 

The (pseudo-)groundstate wavefunction 4>o{x) is determined as a zero mode of A, A(j)Q = 
(I2.14p . The equation reads 



yjv*{x-q)(j> {x - q) = yJv{x + q)Mx + q)- (B- 1 ) 

Let us define T + = \/V{x) e yp ^V*(x) and 71 = \JV*(x) e~ lp \/V(x). Then the Hamiltonian 
(12. 7p is Ti = T + + T_ — V(x) — V*(x). For the QES case, the compensation terms are added. 
It is obvious that the function part — V(x) — V*(x) (plus possible compensation terms) is 
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hermitian by itself. Let us define two analytic functions F(x) and G(x) as follows: 

g*(x)T+f(x) = (/>* (x)Q* (rj*(x)) y/V^)y/V*(x - iy) ~ iy)P(v(x ~ h)) = F(x), (B.2) 
g*(x)T_f(x) = <f)* (x)Q* ( V *(x)) ^VHx)^V(x + iy) <j> Q {x + iy)P(r,(x + iy)) = G(x). (B.3) 

Then we have 

{T + g)*(x)f{x) = <Pl{x + iy)Q*{r 1 *{x + h)) \/V*{x)y/V{x + iy) <f> Q {x)P(r}{x)) = F(x + 27), 

(B.4) 

(T„g)*(x)f(x) = <P* Q (x - iy)Q*(r,\x - iy)) y/V(x)y/V^ - iy) <f> (x)P(v(x)) = G(x - iy). 

(B.5) 



By using (IB.lj) and the 'reality' of c/>o(x), rj(x), P(rf), Q(rj), we obtain 



g*(x)T + f(x) = V(x) Mx) 2 Q{v{x))P{v(x - iy)) = F(x), (B.6) 

g*(x)T_f{x) = V*(x)Mx) 2 Q(v(x))P(v(x + iy)) = G(x), (B.7) 

(T + g)*(x)f(x) = V(x + iy)<j) Q {x + h) 2 Q(v(x + iy))P(v(x)) = F(x + iy), (B.8) 

(T„g)*(x)f(x) = V*(x - iy) o (x - iy) 2 Q{ V (x - iy))P( V (x)) = G(x - iy). (B.9) 

Therefore the necessary and sufficient condition for the hermiticity of the Hamiltonian be- 
comes 

/ (F(x) + G(x))dx= / (F(x + iy) + G(x-iy))dx. (B.IO) 

J Xl J Xl 

Of course it is required that there is no singularity on the integration contours. 

Let C± be the rectangular contours x\ — > x 2 — > x 2 ± iy — > x\ ± iy — > Xx and D-t be 
the regions surrounded by C± including the contours. Under the assumption that F(x) and 
G(x) do not have singularities on C + and C_ respectiveljj^], the residue theorem implies that 
(IB. 101) is rewritten as 
n 

\ (F(x 2 + iy) - F{xi + iy) - G{x 2 - iy) + G(x x - iy))dy 
Jo 

= 2?r n( Res ^(^)- Res xG(x)y (B.11) 

x : pole in D+ x:poleinD_ 

We will mention several sufficient conditions for ( IB. 111) . If F(x) and G(x) are holomorphic 
in D + and D_ respectively, the r.h.s. of ( IB. Ill) vanishes. In the following we assume this. 



2 If there are singularities on the contours — > £2 ± *7 or a;i ± 17 — ► xi, we deform the contours and 
redefine C± and D± in order to avoid singularities on C± . For simplicity we have assumed no singularity on 
C± in the text. 
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case 1 : X\ = — oo, X2 = 00. 

If 4>o(x) is rapidly decreasing (e.g. exponential in r](x)) at x ~ ±00, then (IB. Ill) is satisfied. 
case 2 : x\ = 0, X2 = 00. 

If (f>o(x) is rapidly decreasing (e.g. exponential in r)(x)) at 2; ~ 00, then (IB. lip becomes 
fJ(F(iy) — G(—iy))dy = 0. This is satisfied if F(iy) = G(—iy), y E (0,7). As a sufficient 
condition for F(iy) = G(—iy), we give the following three reflection properties: 



= (f> {x), rj{-x) = rj(x), V*{x) = V(-x). 



(B.12) 



case 3 : X2 = X\ + u (0 < u> < 00). 

The condition (IB. lift becomes (F(xi+u+iy) — F(xi+iy) — G(x\+uj — iy) + G{x\— iy))dy = 
0. This is satisfied if F(x\ + u + iy) = G(x\ +uj — iy) and F(x\ +iy) = G(xi —iy),y<E (0, 7). 
As a sufficient condition for F{x\ + u + iy) = G(x\ + u — iy) and F(x\ + iy) = G{x\ — iy), 
we give the following reflection relations and the periodicity: 

<j>a(—x + xi) = (f) (x + xi), r)(—x + xx) = T](x + xi), V*(x + x x ) = V(—x + x t ), 
<f) {x + 2uj) = cf) {x), r]{x + 2u) = 7](x), V{x + 2u) = V{x). (B.13) 

C Eigenvalues and eigenvectors for an upper triangu- 
lar matrix 



Let A = (ctij)i<i,j< n De an upper triangular matrix, namely a»j = for i > j. Its eigenvalue 
are cuj = an (i — 1, . . . , n). When a^s are mutually distinct, the eigenvector corresponding 
to the eigenvalue Ui is given by the determinant 



ei e 2 e 3 

a>2 — (Xi a 2 3 







where {ej} is the natural basis, (e»)j = 5, 



(C.l) 



■1.1 ■ 
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